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SYNOPSIS 


Continuous litre -reinforced composite materials 
are an important class of modern materials, Tlieir high 
strength to weight and modulus to vreiglit ratios, contro- 
lled anisotropy and formahility make then superior to 
many conventional materials for various applications. 

The strength^ elastic modulus and many other properties 
of such materials are highl37' direction dependent* The 
present studies were carried out to show the feasibility 
of obtaining theoretical optimum parameters of laiilna,ted. 
composite beams and plates subjected to multiple behaviour 
constraints, A number of specific optimization problems 
were solved using the penalty function method of optimiza- 
tion for constrained problems and the Davidon - Iletcher - 
Powell method for solving the sequence of unconstrained, 
minimization problenls. Tests were performed on the 
fib re -reinforced composite plates at different orientation 
angles to find the natural frequencies of vibra.tion for 
simply?- supported boundarj^ edge conditions. 



CHAPTER I 


IHTRODUCTION 

The need for high strength to weight ratio 
materials has increased considerably with the advance 
of technology. Specially in structures where weight 
is a problem, e.g. space-crafts, aircrafts etc., the 
high specific strength materials are desirable. Since 
the directional properties can be controlled very pre- 
cisely in composite materials, they have a higher 
specific strength as compared to the popular isotopic 
materials,' Since the structures are in general sub- 
jected to vibrations, in-plane loading, transverse 
loading etc., the studies have been conducted to find 
an optimum solution with different combinations of the 
above constraints. 

1.1 LITERATURE SURVEY 

The study of optimal design of fibre reinforced 
composites (RRC) has started recently with increased 
application of high performance F.R.C materials for 
efficient design. The minimum wei^t design of lami- 
nates for strength and membrane stiffness was studied. 
extensively by Foye [ 9 ] , Multiple in-plane loading 
conditions were considered and a random search method 
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was used to find the ply orientation angles such that 
the strength and stiffness requirements would he satis- 
fied with the smallest number of plies » Another proce- 
dure for the optimum design of laminates has been 
reported by V/addoups |io]| . Miriimm weight designs are 
obtained considering strength constraints under multiple 
distinct loading conditions. Either the Hill -Tsai or 
the minim-urn strohn criterion of failure nay be used and 
all lamina are ass-umed to behave linearly to failure. 

The search procedure employed is a systematic ’try them 
all' procedure which is quite exhaustive. Both of 
these studies deal directly with discrete number of 
plies and they also treat ply orientations as design 
variables. Waddoups and Co-workers [ill have reported 
a structural synthesis capability for a class of aniso- 
tropic plate structures. Kicher and Chao 02] have 
reported the development of a structural optimisation 
capability for stiffened fibre -composite cylinders. 
Multiple load conditions are considered and strength 
as well as buckling failure modes are guarded against. 
The design variables include stiffener dimensions and 
spacing, fibre volume content and? ply orientation 
angles, Schmit and Earshi have developed a 
laminate optimisation capability in which the thick- 
nesses of material at specified orientation angles 
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are treated as the only design variables . That is to 
say, the materia], to he used would he selected in 
advance from a discrete set of available well-charac- 
terized systems and the orientation angles to he used 
would also he preassigned. By treating the thickness 
of material at each preassigned orientation as a conti- 
nuous design variable the difficulties associated with 
discrete integer design variables are avoided, 

1 ,2 PBESEHT WORK 

In tliis work, the optimum design of symmetric 
laminates is considered by taking the ply thicknesses 
(t^) / the number of plies (N) as design variables so 
as to satisfy certain specified constraints on the 
natural frequencies, buckling load and deflection. 

An experimental study has also been made to 
find the dynamic response of the laminates at different 
orientations since the closed form solutions for different 
orientations are not available , 



CFLA.PTER II 


ANALYSIS OF LAI4INATSD BEA^E AND PLATES 

A laminated ‘beam or plate is composed of an 
arbitrary number of bonded plies, each of which may 
be inhomogeneous having different thickness and ortho- 
tropic elastic properties , As a mathematical approxi- 
mation, the laminate is considered to be an in-plane 
homogeneous and transversely heterogeneous continuum. 

The transverse heterogeneity has a step-wise variation 
in material properties between plies , 

The analysis of deflection of a simply suppor- 
ted beam under uniformly distributed load conditions fs 
done within the frame work of the small deflection theory 
of thin plates and based on the following assumptions : 

1 , The material homogeneity through the plate 
thickness is there, which means that the 
gross laminate properties will apply. 

2. If the Saint Venant's principle is imposed, 
the plate length must be much greater 
than the depth ’2b’ . Also, it will be 
assumed that the plate thickness ’h’ is 
unity, much smaller than '2b’ . 
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Thus the analysis is, theoretically, appli- 
cahle only for a long, tliin, homogeneous anisotropic 
plate . 

The analysis of dynamic behaviour of the 
laminated pla,tes is also done within the frame irork of 
the small deflection theory of tMn plates and "based 
on the follovri-ng assumptions ; 

1', Each layer is homogeneous, linearly elastic 
and under a state of plane stress. 

2. The plate thickness 'h' is small in relation 
to the lateral dimensions, and the transverse- 
deflection w is small compared to h. 

3« Plane, normal cross-section "before deformation 
remains plane and normal to the reference 
surface of the plate during deformation; thus 
transverse-shear deformation can "be neglected. 

Under the above assumptions, the problem of the trans- 
verse vibration of the laminated plates without in- 
plane forces becomes that of' a homogeneous plate once 
a constitutive equation for the laminated plate is 
obtained , 
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2.1 LAMINATED PLAIES 

In this section we descrihe the general notations 
and constitutive equations for the laminated plates with 
appropriate simplifications for using in further analysis. 

2.1.1 General Laminated Elates 

A laminated plate consists of 2W constituent 
plies of different orthotropic material and thickness 
with arbitrary ply orientations . The structure axes of 
the laminate are denoted by x, y, and z, and for each 
ply the t\<ro symmetric axes of material in the plane of 
the ply by L and T. The ply orientation angle 0 is 
the angle between the x - and L - axes, and measured 
in the counter-clock -wise direction from, the positive 
X - axis as shown in figure 1 . The structure coordi- 
nate system is located so that the X - y plane (z = 0) 
coincides with the middle surface of the laminate . 

The plies are numbered from 1 to 2N from the bottom 
ply to the top ply with the bottom surface of the ith 
ply located at the distance from the middle surface 
as shown in -figure 2. The ply orientation of the 
laminate will be denoted by (9.|, 02 ••• ). 

For orthotropic plies under a state of plane 
stress, four independent elastic constants are requi- 
red to completely characterize the elastic behaviour. 
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These principal elastic constants are the Yonng’s 
moduli in the two principal directions and shear 
niodii].us 5 GpT? and the major Poisson’s ratio, 

The constitutive equation, which characterizes the 
structural hehaviour of the laminate, can he obtained 
once the principal elastic constants, thickness and 
the ply orientation of each unit ply are known. 


The constitutive equation is the relations 
between stress resultants {n} and stress couples 
{m} and the reference surface strains ^^S^jand 
curvatures {OCV , By defining the stress resultants 
and stress couples according to the linear theory of 



the constitutive equation for the laminate can be written 
in the following matrix fom : 



^12 

^13 I 

’^I'l 

^22 

■*•23 I 

®12 

A23 

A ^ 

_ 33 _ j 

13 

^ 12 "' 

®13 i 

ll 

®22 

®23 i 


®23 

Bqo ! 

I 3 


®12 ^13] ^oz 

®22 ®23 ^oy 

®23__ ®33 ^ ^oxy^ 

^12 ^13 

^22 ^23 

^23 ^33! ^'^xy 

-j ^ ^ 


( 1 ) 
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where, [a] is the in- plane membrane stiffness 
matrix, lb ./in. 

L^J 'tlie stiffness coupling matrix, lb. 

r n 

and ^DJ is the flexural stiffness matrix, lb. -in. 

Since the reference surface (z = 0) is talcen to coincide 
with the middle surface, the stress couples are computed 
with regard to the middle surface, and the reference 
surface strains are the middle surface strains. The 
reference surface could be chosen to be other than the 

f *~T 

middle surface. Consideration of gBJ matrix in the cons- 
titutive equation indicates the coupling phenomena bet- 
ween stress resultants and the curvatures and between 
the stress couples and the middle surface strains which 
do not occur in the theory of homogeneous plates. 

It is convenient to express M and in terms 
of X and N, for the bending ansilysis of laminated 
plates. Solving the first three equations in (1) for 
in terms of X and N, and introducing the results 
into the remaining three equations yields , 


= [a*3 n + [b*] 

M = [c*j K + [d*3 

where = [aJ 

[b*] = - [a] [B] 

[c*3 = Fb] |a] '1 - 

.d [D*] = [D] -■ [b][a]-^ [B] 
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It is noted that while 1_A | and [_D’’ | are 

symmetric matrices, the same need not he true for |b''J 
and [p J , If the laminated plate is free from the 
stress resultants (no external in-plane force), the 
stress couples and the middle surface strains become, 


[b*] '< 


r sK*] 

The matrix |_D j is called the reduced stiffness matrix. 


2.1.2 Special Laminated Plates 


Por the laminated plate having an elastic 
symmetry about' the middle surface, the coupling stiffness 
matrix is 'a null matrix, and the bending a,nd membrane 
behaviours are uncoupled. The reduced stiffness matrix 
is, therefore identical to the flexural stiffness matrix, 
and the middle surface strains due to bending vanish. 


For the purpose Of this study, the lamj.nate is considered 
to be orthotropic if the laminate is uncoupled and its 
membrane and flexural stiffness matrices are orthotropic, 


i .e . -1 


^13 " ^23 


0. An orthotropic laminate 


has three planes of elastic symmetry:’ the middle surface 


and the two planes whose normals are along the x and y 


axes 
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Por coupled laminates, the sufficient condi- 
tion that the reduced stiffness matrix is orthotropic 
= ^22 = 0) is 


(a) 

Ai 3 

- A23 - 

0 



D ^3 

= ^23 = 

0 



®13 

= B23: = 

0 


or 





(b) 

^13 

= A23 

0 



^13 

= D23 = 

0 



®11 

= ®12 = 

^22 “ ^33 ~ ^ 


The elements 

® 23 ’ ® 31 * 

B32 vanish also for 

the 

condition 

(a) , 





For the 

laminated plate consisting of 

const! 


tuent plies of the same orthotropic material and 
thickness such as a laminated fibrous composite plate, 
the constituent equation becomes only function of the 
ply orientations once the principal elastic constants 
and thicloiess of ply materials are specified. Depending 
upon the manner in which the individual plies are 
oriented X'd-th respect to each other and arranged about 
the middle surface of the laminate, some of the elements 
of the membrane, coupling and flexural stiffness 
matrices vanish as follows : 
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1 . 


2 . 


3. 


If the plies are oriented sjnnme trie ally about 
the middle surface, (9^, ^3? ®2’ 

All E - , = 0 


If the plates are oriented anti -syinme trie ally 
about the middle surface, (9^, 9*2, 9^, " ^35 
- 92, - Q-j), then 


■^13 ^ ‘^23 ^ ° 

^13 "" ^233= 0 
®11 ~ ^12 ^ ^22 = ^33 " ^ 


If some of the plies are oriented symmetrically 
and the others are oriented anti -symmetric ally 
about the middle surface (9^, 92, S3, -■ 02> 6-])) 
then 


if. 


B 


11 


= r^12 = ^ 


22 


= B 


33 


0 


If some of the plies are oriented symmetrically 
about the middle surface and the others are 
oriented either at 0° or 90° (9^, 0*^, 9^, 90°, 9^), 
then 


5 . 


®13 " ®23 " 


0 


If the plies are oriented at either 0° or’’. 90° 
regardless of symmetr3r about the middle 
surface, then 


A 


D 


B 


13 

13 

13 





0 

0 

0 
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Wlien the laciinate has isotropic plies, such as glue 
layers in fiherous composite plates, in addition to 
the orthotropic plies as* its constituents, case {^) 
is uneffected, and case (l) to case (^) are also un- 
effected provided the isotropic plies possess identical 
thickness and elastic properties, and are located 
symmetrically about the middle surface , 

2,2 ANALYSIS OF LAMNATED EEAI-S 

In this section ve derive the deflection compo- 
nents u and V in X - and y - directions for a simply 
supported beam under uniformly distributed load (see 
figure 3) • The approach is same as the one advanced 
by Hashin [l^J Tor anisotropic beams in a state of 
plane stress, 

2,2.1 Boundary Conditions 

The boundary conditions to be satisfied by the 
stresses and loading are as folloi^rs : 


<> 

(x, 

-b) 

= 0^ (x, -b) = 0 


(z, 

b) 

= 0 


(Zl 

b) 

= q. 

M. 

o 


0 


M. 


0 



(3a) 



13 

= 0 
= . 0 
= ~ql/2 

(3Td) 

where cr and o~ are normal stresses in the x - and y - 
directions respectively, ^■>ry is the shear stress in the 
xy plane, P and II represent the force and the moment, 
respectively, 1 and 2 representing x - and y - directions 
and o and 1 representing values at origin and at length 1* 


P 


1o 


11 


■2o 


21 


2.2,2 Stress Function 

In the absence of body forces, the plane equations 
of eqirilibrium will be Identically satisfied by the Airy 
stress function U (x, y) defined by 


_^2 


cr 

X 

6 u 


“ 2 



o y 


""y 

3^ IJ 

- - 2 
d X 

C^) 


2 


t. 

d U 


'-tcy 



3x By 



and the strains £ . . are given by 

J- J 


where the compliance matrix, is the inverse 

of the stiffness matrix, K- when S.,.. is reduced 
to variables mth two subscripts as shown on page 15 . 



The only survi'ving compatibility eq-uation is 


3 ^ f 


= 2 


dx 5 y 


substituting equations ( 4 -) and ( 5 ) into equation (6), the 
differential equation for U can be obtained as 


^2222 ^’xxxx *" ^ 2 ^^iipp ^ ^1P1P^ 


2212 ’xxxy 


1122 


1212-' "»xxyy 


^1112 '^’xyyy ^1111 "^’yyyy 


where partial differentiation is denoted by a comma, 


Using equations ( 3 (a)) and ( 3 (b)) the boundary condi- 
tions for the stress function can be stated as (13I 


U (x, -b) = 


(x - \) o: ( y 5 - b) d^ =0 8(a) 

/ y > i 


U (x, -b) 


( 5 . - = 0 


8(b) 


U (x, -b) 


t^(y, -b) dt . 0 


8(c) 


u (x, b) 


(x -> ) o- ( y , - b) d^ 

J ^ / i 


-2 b 




+ P23_ (1 - x) + b 


4. ( 5 - x) <r ( ^ , b) d; 


■x 


U (x, b) = ^ (21 - b - 3 x - xVl) 


8 (d) 



u (x, b),y 


~ ” I ^ ”• "b) 

n 


+ p 


11 



■ 15 


.8(e) 


In tlie following discussions, to avoid complicacy, 
the following notation will he used for. the compliance 
matrix elements % 


^1111 ®11 ’ ^2222 " ^22 » ®1112 ~ ^13 
^1122 ~ ®12 > ^1212 ^ ^33 


The other elements can he defined in a similar manner. 

Let us assume the stress -function to he defined 
as follows ; 


U (x, y) 



m=0 


^B-3 

IZ 

n=0 




such that m + n M + 3 and M is the highest power of 
X in equations (8 (a} to (8 (e)) 


. M = 2 and m + n 2 + 3 = 5 * 

.. U(x, y) = ^00 + Cq^ y + ^Q2 ^ ^Q3 ^ ^h':- ^ 

+ Gq^ y^ + C^q X + xy + C^2 

+ xy^ + xy^ + G^q x^ 

+ y?J + 0^2 ^ 7 ^ + C22 (9) 
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Using compatiUllity equationSj stress-strain 

relations and stres s -function , we get the recursion 

relations for G as 
mn 

^22 (m + 1) m (m - 1) ~ ^ ^23 ^ ^ 

* (m - 1) (n - 1) + 2 + 2 m (m - 1) n 

^mn - ^ '^3 - 1) (n ^ 1) n (n - 1) 

+ (n + 2)(n + 1 ) n (n - 1 ) C^_2^n+2 ~ ° 

when m > 2 and n > 2 

The total number of coefficients is given by 

IQH 

S = 1 Cm + 1) (M + 8) = 1 X 3 X 10 = 15 
and total number of recursion relations is given by 

R = I M (H + 1) = 1 X 2 (2 + 1) = 3 

t 

. . For m = 2 , n = 2 we have from equation (10) 

(S12 + 2 0^2 + 3 Gq^ - 3 0^2 = ^ 

For m = 2, n = 3 "we have from equation (IO) 

(S12 + 2 ^ 23 ^ *^23 "* ^ ^13 ^l 4 ^ ^11 ^05 ~ ^ (12) 

For m = 3j n = 2 we have from equation (10) 


^13 ^23 ^11 


0 


( 13 ) 
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By evaluating the partial derivative of U, (x, y) 
id. til respect to y at y = b and equating coefficients 

of like powers of x in this equation and equation 
^8(e))we get 

^01 *^02 *^03 ^ '^Qif + 5 Cq^ b "^ = - qV 2 

(14) 

°11 213 C.|2 + C.|2 + C-llj. = 0 (15) 

Cg., + 2b 0^2 + 3b2 O23 = 0 (16) 

Again evaluating the partial derivative of U (x, y) with 
respect to y at y = -b and equating coefficients of 
like powers of x in this equation and equation (8( c ) ) 
we get 


°01 - 2^= S 2 °03 - d Cgj =0 (17) 

°11 ■ 2b + 3b2 C.,3 - Vb^ = 0 (18) 

- 2b 0^2 + 34^ C22 = 0 (-19) 


Now evaluating U (x, y) at y = b and equating- 
coefficients of like powers of x idth those in equation 
(8(d)).^ %fe get 


°00 °01 * + b3 + b^b + b? = ■■!^s(21-b) 


02 


03 ^ ^04 


"05 “ 2^ 
( 20 ) 


^10 ^ ^11 ^2 ^13 ^ 


^4 = - 3ql./2 

(21) 

20 ' ^21 ^ ~ q/^ (22) 

Similarly evaluating U (x, y) at y = -b and 


C20 "*■4 02.^ + b C22 t)^ C22 — ~ q/2 
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equating coefficients of like powers of x ’td.tli those in 
equation (8(a)), we get 

^00 “ ^ ^01 ^02 - ^03 ^05 = ° ^ 23 ) 

^10 " ^11 * ^12 “ S3 ^14 " 0 ^ 24 ) 

C 20 - + h^ 0^2 - 43 C 23 = 0 (25) 

A simultaneous solution of equations (1 1 ) to (25) gives 
all the 15 coefficients, to completely define the stress- 


function U (x, y) . The coefficients are given helow; 

Cqq = ql ((41 - b) + 3 b) / 8 (26) 

Gqi = - ql 01 + b)Ah + qb (8 - (S^^ 333 ) 3^^)740 3^^ (27) 

3q 2 = ■* ql/8b - 33^3ql/4S^^b ( 28 ) 

C 03 = ql^ / 4b3 - q (8 s 2 ^ - (3^2 + 2 3^^) 3^^) / 20 b 3^^ ( 29 ) 

*2q 1^ — 3 3^2'.ql / 8 3^^ b^ (30) 

S5 “ q ^13 “ ^®-]2 '*■ 2 ^33) S-Il) / ^0 32,^ b^ (31) 

C-IO = " 3 ql / 4 - 3^^ qb / 4 3^^ ( 32 ) 

I'll = 9 ql / 8 b (33) 

C^2 - “ fl / ^ (3^) 

0^3 = 3 ql / 8 t3 ( 35 ) 

C'li,. = S ^3 q / 4 3^^ b^ ( 36 ) 

Cpo = -q / ^ (37) 

^ 2-1 = - 3q / 8 b ( 38 ) 

0^2 =0 (39) 

C 23 = q / 8 b^ • (40) 
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2 . 2*3 Deflection Eque.tions and the Stress Components 
The deflections are given hy 

(^!-1 ) 

and V = S^2 ^’y ^ ’^22 •' ^ ^23 U,^ f^Cx) 

(^-2) 

The arhitro.rj^ functions f^ (y) and f^ (x) are given by 
using 

“^'xy “ 2 xy ~ 2 ^'^>y ''>x^ ~ ®13 '^’yy ^ ®23 ^’xx 




( 43 ) 


and three boundarj'- conditions , namely u = 0, v = 0 at 
X = 0, y = 0 and v = 0atx = l, y = 0. 


Finally, we get the deflection equations as 
u (x, y) = 2 (S^^ Cq2 - x 

+ Gi- (S^2 Cq2 - + S 22 G^q) 

(3 C^2 + ^^12 ^ ^33^ ^21^ ^ • 

~ _(C^2 + ^23 ^ *^^11 S2 “ 2 ®13 ^21^^^ 

+ 2 (3 Cq2 - 2 S^2 S2 ®12 ^21^ ^ 

^®12 S2 + 2 (3 Gq2 - <^S^2 2 C^2 

+ 2 S23 C^^)) y2 + 3 G^2 

4- 6 (2 S^^ xy^ + (8 S^^ Cq^ - ^12^13 

7 ^ ^33 ^13^ 7^ + 2 S^^ C23 y +6(S^^ 

S^3 G23) x^ 7 ^ + 2 (10 Cq^ - If 

+ C^^j) xy^ (l+if) 
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V (x, y) = ((3 - h + if S,,) C^^) 1 


13 12 


^3^ 21 


+ (C^^ + 0^3 1) 1^) X 

+ 2 (S^2 ^02 " ^23 ^11 ^22 ^20 ^ 

- (3 Cq 3 + (S ^2 + ^ S 33 ) ^21 ^ ^ ^13 ^2^ 

+ 2 (S^2 C']2 - 2 3^3 C^^) xy + (3 3^2 ^o3 

- 2 8^3 C^2 + ^22 C2^) y^ - C^3 x^ 

3 ^•]2 *^13 ^ +2 (2 S^2 3q]_^ - 323 ^13^ 

- ®11 ‘=23 3 S^2 C23 x2 y2 


+ 4 (S.o - S , C„,) xy^ + (5 .s.„ C, 


12 14 23 23 “ 

” ^23 ^14 ^ 2 ^22 ^23^ 


12 "05 


(45) 


The mid-plane displacements are given by 


u (x, 0) = 


2 Cq2 - 3^3 + S^2 ^ 20) ^- 


+ C^2 - 2 3^3 C2^ ) x'^ 


V (x, 0) = 


(3 3^-j Cq3 - 4 S^3 0^2 + 1 2-11 35^3 


4 (8^2'*“^ ®33^ ^21 ^ ^11 

- (3 3^^ Cq3 - 4 S^3 C^2 ■*■ (3^2 ^ ^33^^21^ 


S.. C.2 -■ x^ 


’11 ^^13 


'n ^23 


(46) 
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The stress components for a simply supported iDeain 
under uniforaily distributed load as given in reference |^7j 


are (after transforming to the 
in figure 3 ) : 


II 

(k) 

A 

B M 

— X 

(x) 

II 

{10 

h2 

A 

B M 

— X 

(x) 

and 

^(k) 

(lO 




xy : 

= 

A 

E M 

“ X 

(x) 

where 

A = 

All 



1 

C\) 

"'Hi 

hi 


B = 

-li + 
hi 

y j 




coordinate system shovm 

(if7) 

(^7(a)) 

(2+7 (b)') 


K. . represents the stiffness at any angle 9 and the 
1 3 

superscript k indicates the ply* 
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2.3 ANALYSIS OF LAMINATED PLATES 

In this section we derive the expressions for 
the transverse deflection, the critical huciaing load 
and the frequency of natural vibration for a simply 
supported laminated orthotropic plate within the frame 
work of the assumptions made in section 2.1, and at 
the beginning of this chapter. 


2.3.1 Transverse Deflection Analysis of 
Laminated Plates 

The classical laminated plate theory |j , 2, 3^ 
reveals a coupling phenomenon between bending and middle 
plane extensions which was not found in the linear 
theory of homogeneous plates . This phenomenon occurs 
through the constitutive relations given in equation (1). 
Due to limited existence of solutions to boundary value 
problems, the effect of coupling on the mechanical 
behaviour of composites is not well understood. In this 
section the effect of coupling on the response of simply 
supported rectangular laminates of special construction 
under transverse load is studied | . 


Closed form solutions can be obtained for 
coupled laminated plates in which j^Aj and ^d]| 
matrices of equation (1) are orthotropic (i.e. A^^ = 


•^23 ~ ^13 


and either ^22 ®13 ®23 



23 


are the only non-zero elements of the [b 1 matriii. 

Using a standard x, jj z coordinate system, the govern- 
ing equations for unhalanced cross -ply and angle ply 
plates in terms of displacements can he expressed as [^3 J • 


A 


11 


u, 


'XX 


+ A 


u, 


33 ‘"’yy 


4* (A 


12 




A33) 


V 


J-'r 


xy 


- B 


11 




XXX 


- 3 E 


w, 


- B, 


'13 '"’xxy '■ ^■"23 "'^’yyy 


= 0 


(^8) 


(A.„ + A^-,) u' 


^12 


33' 


'xy 


4- A 


V 


+ A 


V 


- B, 


33 ''»XX ' ^"22 ''’yy ^22 ^'^’yyy 


- E 


13 




XXX 


-^^23^'^’yyy = ° 


(^9) 


^11 ^'^’XXXJC ^ ^ w,^yy. + D 22 ^hjyjj 


’■^11 ^’xxx “ ^22 ^’yyy 
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XXX 


’xxy' 


- ^23 '3 




(50) 


•where u°, v° and w denote the components of displacement 
of the middle plane along the x, y and z axes, respec- 
tively, a camma indicates partial differentiation. 


Consider a rectangular cross-ply composite 
of dimensions a, h having an even number of plies with 
each ply ha''/ing the same thickness alternately oriented 
at 0*^ and 90° to the x - y axes of the plate. . 
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^22 ~ ^22 = ^11 > - B 


(E - E„) ^ 

— -4» ■ 1 -1^ 7T 

22 8 Em W ^ ^11? 


^13 ~ ®23 - ° 


(51) 


Assume that the plate is simply-supported in such a 
manner that normal displacements are adriiissible, but 
lateral (tangential) displacements are not. Thus the 
foil owing boundary conditions are appropriate : 

"w (0, y) = (a^ y) = y Q) _ ^ 

= (0, y) = (a, y) = (x, o) - 

= My. (x, b) = 0 ( 

v° (0, y) = v° (a, y) = u® (x, 0) = (x, b) 

= (0, y) = (a, y) = (x, O) 

= Ny (x, b) = 0 ( 


A solution can be ootained by expanding the 
transverse load q into a double Fourier series as 


II 

OO 

^mn 

Sin 

m TI X . n K. V 
— r — sin 

a -u 

(54) 

m=1 

n=1 



u 


and solutions of 

the 

form 




II 

o 

£ 

A 

mn 

CDS. - Sin —Tv X 

a h 

(55) 

m-1 

n=1 





0 & 

^ = E 

00 

E 

B 

mn 

sin 

mK X , . , nil V 
,a -b' 

( 58 ) 

m=1 

n=1 
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CD 

£ 

m=:i 


n=l 


mn 


sin 


mTt X 


Sin 


nTX V 


(^) 


The boundary conditions, equations (52), (53) are iden- 
tically satisfied by the displacement field, equations 
(55) to (57) • Substituting equations (5^) to (57) into 
the governing equations (48) to (50), taking equation 
(51) into accour.t, and solving the resulting simultaneous 
equations, for the Fourier coefficients, yields- 


mn 


‘^mn ^11 ^ 4 . 2 2 t 52 

' n m + j m n R 

-pf 3 p 33 n 

mn 


+ (A^2 


(58) 


B 


Imn 


mn 


k3 


D. 


11 ’^'- " lu 

((^.-^2 -^ 33 ^ “ 


mn 


+ A^^ n^ + A^^ P^^) 


(59) 


mn 


n 

%n ,,2 , fl „2 p2^ 

m ^ "^33 ^ It y 




mn 


(A^^ m^ -^11 "■ ^^12 "^33^^ 


( 60 ) 


VJhere R is the length to xd.dth ratio of the plate, a/b. 
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and D 


mn 


( ( (A^ ^ '^■^'11 

- (A^2 + ^33)^ (m^ + n^ R^) 


+ 2 (D^2 + 2 ^33) “1^ R^) - Rii^-^.'l'l R 


2 2 „2 


* (m^ + n^ R^) + 2 (A^p + A 33 ) n^ R^ 


+ -^33 ^ ^ ■*■ R^ ) ) ) 


( 61 ) 


If the load q is a sinusoidal distributed load 
as given by 


q (x, y) = 


q^ sin 
^o 


TX X 


sin 


ILz 

b 


( 62 ) 


and the plate is a square lanrLnated orthotropic plate i.e. 
coupling matrix [bJ is identically zero, then for a unit 
load q^ i.e. q^^ = 1 the maximum transverse deflection at 
the middle of the plate is given by 


4 




max 


2 + D ^2 + 2 


(63) 


because m = n = 1 in equation (62). 

2 . 3.2 Buckling Analysis of Laminated Plates 

The differential equation describing the bend~ 

/ 

ing of an orthotropic laminated plate subjected to an 
in -plane edge loading (neglecting any bending-membrane 
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coupling term) is f6, 7 


^ + 4 D 

11 ’xxxx 


-i:? + 2 w,. + ^ D„-, w. 

13 ’xxxy 3 23 ’xyyy 


+ Doq ¥, + N ¥« + 2 I'J w« + w 

22 ’yyyy x ^ ^xy ''"’xy y 


'yy 


0 


(64) 

For the orthotropic case with the principal axes 

parallel to the sides of the plate, and with ?I = N = 0, 

’ xy y ’ 

equation (64) hecomes 


^11 ^^xxxx ^ ^3 ^22 ^'’^’yyyy ^ \ ^»xx ~ 


(65) 


where + 2 

3 12 -33 


In the case of a rectangular plate with all the 
four edges simplyrsupportedj the mode shape is assumed 
to be 


w 


A 


mn 


sin 


mKx 

a 


sin 


nlly 

b 


( 66 ) 


substituting equation (66) into equation (65), the critical 
buckling load can be found to be 


(N ) 

x'^cr 


TC 


b" 


D 


m^ b^ 


L 


11 


+ 2 F 


2 ^ T^ 

3 n + Dgj 




( 67 ) 


where m and n are chosen to make (Nj^-);.^ a minimum. 
Since n appears only in the numerator, the minimum 


occurs at n = 1 
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For an isotropic plate 
and equation 67 lecomos 


D 


22 


= D 



E = D 

12(1 


^Vcr 


D , ml a N 2 

■^2 a mb 


( 68 ) 


and, for a square isotropic plate. 

Of ) = i-iiLD 

X cr ,2 


( 69 ) 


In case the plate has unbala^iced configuration 
the stiffness coupling matrix [b") \^ill he non-zero. If 
the coupling matrix has non-zero terns, the preceding eigen 
value formulation can still he used if the effect of the 
coupling terms can he approximated hy adjusting the flexu- 
ral stiffnesses. This can he accomplished hy first consi- 
dering the force -deformation equations (1). 

These equations (1) can he solved for -[ll) hy 
solving the first equation for and then substituting 

its va3.ue into the second equation as follows j 


W = [A],{£j - [B]lX} 

{£0} = W' {”} * W'’' 

{m} = [B] (£„} - [D] 

{m}= [b] [a]--' {tr}. [Biyh Mh} 

{11} = ( [D] - [b] [Af Vb] (-X) +[BjrAl''{n} (70) 
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(71) 


{M}'= [d*J{-X} . [BjrA]-^ U} 

where 

[D*j = [d] - [b] [a]-'’ [ BJ (72) 

The reduced stiffness matrix, can easily 

he calculated and used in place of TdJ in the buckling 
predictions . 

The effect of the menihra-ne -bending coupling terms 
on the buckling load has largely been neglected in the past, 
VJhile it is quite clear from, equation (71) that coupling 
will have very little effect in most cases because of the 
relative magnitudes of the[_Al and [bJ matrices, it has 
been found, in some cases, the coupling terms cannot be 
ignored. 


2.3.3 . Vibration Analysis of Laminated Hates 


The follomng discussion on the transverse 
vibration without in-plane forces are restricted to the 
rectangular, orthotropic laminates and coupled laminates 


with D 
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= 0 for which the theory of homogeneous 


ortho tropic plates is applicable. The effect of material 
damping, the energy-dissipation properties of a material, 


is not included ip. the analysis ♦ The edges of the plate 
having sides of length a and b are assumed to be simply 
supported. For most cases in the vibration of rectangular 
plates, the modal lines of the normal mode are parallel 
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to the sides of the plates, and the normal mode is there- 
fore defined hy m and n, the umber of nodal line in 
the X and y directions, respectiveljr. This number 
does not include the nodal line at the plate edge, i.e, 
the simply supported edge. 


The differential equation of motion for the 


plate element subjected to internally Induced moments and 
externally applied force P as shoun in Fig. 5? neglecting 
body force and rotary inertia, is 



^2 . 2 

o M ^ M 

xy ^ y 

dx 3y 9y^ 


p h 



(73) 


where ^ is the average mass density of the laminate. 
Since shear deformation is neglected, curvatures can be 
expressed in terms of deflection as 






3^ ¥ y _ o 3 w 
3^2 » ^xy - ^ ■^To^ 


Therefore, the moments for the orthotropic 

laminate given in equation (1) in expanded form become, 

2 . 2 




M 


(j) ¥ + I) ^ 

^^11 3^2 ^^12 3^2 


) 


- (D 


M 


xy 


■n2 

O 0 W 

■^33 Sy * 


2 .2 

+ ^ ) 
y 


3 W -r^ 3 W 

12 3^2 22 3„2 


Introducing the above moment expressions in equation (73) 
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yields the following differential equation of motion for 
the orthotropic laminate in terms of deflection; 




'h - 5 + 

d w , j. a \T 

'dx^ 22 ^ 


jj ■'■,2 

V , Cf ^ 
p h 


at-^ 


= p 


(7^) 


The boundary conditions to be satisfied, for the simply 
supported edge whose normal direction is y, are 

0 


w 


M 


7 


2 2 

0 or B ^ = 0 


( 75 ) 


For the coupled laminate, the moment expressions given in 
equation (2) lead to the above equation of motion and 
boundary conditions in which D. . are replaced by the reduced 

stiffness D. .. 

13 

Tlie differential equation of motion and the 
boundary conditions given in equation (7^) and (75) are 
identical to those of a homogeneous orthotropic plate 
except that the coefficients appearing in equation (7^) 
and (75) are the average mass density and the flexural 
stiffnesses (or the reduced stiffness) of the laminate. 

Since the coefficients do not affect the form of the 
solutions, the solutions available for the homogeneous 
plates are applicable for the laminated plate if the 
average mass density and the proper stiffnesses (Dj_j for 
the orthotropic laminate and D. . for the coupled laminate) 

IL 3 
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are used. In the following discussion, unstarred D. . 
will be used for simnlicity. 


The solution for the free vibration of the 
plate which is described by the homogeneous form of 
equation (7^0 leads to normal modes and corresponding 
natural frequencies. No difficvilty is encountered in 
effecting a separation of spatial and time -dependent 
variable for the homogeneous form of equation (7^) . 


Setting w 




m 


n 


mn 


developes that is a liaiTiionic function of time with 
circular frequency P and that VJ must satisfy the 
equation 


'-f 


3 W 


D 


mn 


d VJ 




+ 2 D 


mn 


if 

3 W 


3 h2 


mn _ 




oj 


TT - P ^ Pi: 


w„ 


mn mn 


( 76 ) 


A further separation of as a product of a function 
of X and a function of y is possible for plates having 
two opposite edges sjjnply supported. 


As suming 




= sin 


(m ■» 1 )'n X Y 


’mn a mn 

which satisfies the boundaiy conditions at edges x = 0 
and X = a, and substituting this into equation (76) 
yields the following ordinary differential equation 


for Y 


mn 
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X. 


D 


22 




mn 


(m + 1) d^ X 
p D - Kin 

3 .2 


a 


^ / vli- II 

, i (m + 1) ■ ' 

^ L ^ 


a 


dy‘ 


h p2.„„ i X. 


-] j — -^ rnn i mn 


0 


The general solution of equation (7) is 


(77) 


^im = Sm a 


^mn y ■ y 

— + C sinh 

^2mn a 


+ c, ^ cos "heX, + 

3mn a 


Sin 


'i : ^ y 
rmn ^ 


hrnn 


(78) 


■where 


. n* In -n .^2, 1/2, 1/2 

0 = ^ b _ D D + ^ 

< ft 11 12 + 


mn 


'll)- ■ 

'■ mu 


yD 


22 


(m + 


iiiJ.x :. -r\ ‘ 

xxy“i 3j 


(m ->- Dy . j 

7^22 


❖ ifP-, 1/2 

r f r P h D^p a p i 1 


( 78 a) 


The parameters and are related by the equation 


mn 


D. 


0^ - =:2(m'+’l)^TT^ t?3— 

^mn rmn <^'vm i;h 
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The constants are determined by four 

homogeneous boundary conditions, two at the edge y = 0 
and two at the edge y = b, given in equation (75) • i’or 
homogeneous boundary conditions, only the ratios of the 
^iran determined through the use of any three equa-^ 
tions, while the fourth boundary condition leads to the 
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frequency equation. The roots of this equation are the 
eigen values uhich tiie solution is valid, 

lor each thei"e is corresponding eigenfunction which 
constitutes a normal mode of vihration of the X-^late. The 
complete solution is obtained by supeirposition of all the 
normal modes ; i . e . , 


w = n S sin - ( - a -t -O - ix ( , 

m n a mn ^mn 


n, m = 0, 1 , 2 


( 79 ) 


The circular frequency p^^^^ is a natural frequency of the 
plate and relates to the roots 0^^ of the frequency equa- 
tion by equation (78 (a)) and is given by 

P h 2_2 h if ^ q 1/2 

_ L h? ^ ^ ^ °11 J 

^mn „2 /— 5 

Vp h 

j 

The fre ;^uency equation for a plate simply suppor- 
ted on all the four edges is given by 


sin TI) R 
r mn 


0 , 


ih - 

^mn R 


the roots of which are 
where R is the aspect -ratio of the plate. 


(n+l)'7(‘ a 
b 


K 


■mn 


Therefore the natural frequencies are given by 

2 r- /i- ..0..0 ..vif 1 

1 

22 J 




hi * 2 

w*. 9/ 


(m+1)^ (n+1)^- 


2 ^2 
a b 


D. + D, 


1 / 


( 80 ) 
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and by 


mn 


sin 


(n^1)K 

Id 


Z 


(81) 


The nodal patterns of a normal mode consists, 
of straight lines parallel to the plate boundaries of 

( a) ni lines given "by sin = o in the 

a. 

X - direction. 

( ^ 1 I ' i 

(h) n lines given hy sin ^ = 0 in the 

y - direction. 

The.i'^ave forms in the x - and y - directions 
for SSSS horuidary conditions are sine waves and are 
independent of the values of the elastic moduli. 



CHAPTER III 


FORl'fJLATIOlI AHD SOLUTION OF OPTIMSATION PROBLEM 

When a means for predicting the behavior of any 
design vnthin a particular design concept is available, 
limitations on the performance and other external cons~ 
traints on the design can bo stated, and an acceptance 
criterion can be established, it is possible to cast the 
design modification problem in the form of a mathematical 
programming problem, -wiiich can be stated as follows s 

Find X which minimizes F ( X ) 
subject to the constraints 

gj- ( X ) 4 0, j = 1, 2, , p (82) 

where X is called the design vector consisting of the 
design variables x^, X2j ... x , F (X) is called the 
objective function and g. ( X ) is the j constraint 

J 

function. 


3 . 1 FROFLU:: state: ENT 

3.1.1 For Laminated Beams 

In this case, the design determination of the 
cross-section of a simply supported beam under uniformly 
distributed load is considered for minimum volume when 
the load and the elastic properties of the beam are 
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specified. CoxivStraints are specified on ilie noimal and 
shear strengths and on the maximtmi deflection. 

More specifically, the following optimisation 
problems were considered. 



3. Upper limits on the normal and shear stresses in 
laminate as X, Y and S respectively. 

Load distribution on the beam. 


5 . Span ’1' of the beam, width of the beam h = 1 . 
Find : A set of i = 1, 2, ..., U, such that 


1 . 

Volume 

= 21 23 V . — minimum 

i=1 

( 83 ) 

2. 

cr 

X 

[ max 

4 X 

(8lf) 

3^ 

Oj 

i ^ 

: mas 

< I 

(85) 

h. 

4 S 

xy max ^ 

(86) 


5 . 0 for i = 1, 2, ..., U ( 87 ) 

+ If, as is often tlie case", the laminate is to be made from 
a single ts^pe of fibre composite material, only one set 
of material properties need be specified. 



It should he noted that ec[iiation (83) states the minimum 
volume objective, equations (84) to (86), implement the 
stress constraints and equation ( 87 ) requires that the 
thicknesses be non-negative. 

Problem II ; 

Given ; 

1 . The available orientation angles ' 0^^ and the 
corresponding material properties , 

and for i = 1, 2, N. 

2. N the number of available orientation angles. 

3 . Upper limits on the normal and shear stresses in 
laminate as X, Y and S respectively, 

4. Upper limit: on the maximum transverse deflection 
of the beam. 

5 . Load distribution on the beam. 

6. Span ‘1’ of the beam, width of the beam h = 1, 


Find; 

A set 

of i = 1, 2, U, such that 

N 

n 

(88) 

1 . Volume ; 

= 21 /Lj 4 . — minimum 
i=1 ^ 

2 . 0-^1 

max 

4 X 

( 89 ) 

3» cr j 


X 

( 90 ) 

y 

max 
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5 * ^ V ' ' (92) 

max ^ o 

6. 0 for i = 1, 2, ..., IJ (93) 

It shoifLcl te noted that equation ' (38)- states the minimum 
weight objective, equation (89) to (9I) ^implement' the 
stress constraints, equation (92) implements deflection 
constraint and equation (93) requires that the thickness 
he non-negative. ' ' ■ 

It can he noticed that the ahove two problems 
are similar to the one stated in equation (82). 

3..1:,2 For Laminated. Plates 

In this work, the ply thickness (t^) / half the 
number of plies (h) are treated as design variables for 
minimum weiglit design of laminates for preassigned values 
of frequency, buckling load and- ’.deflect! on. ■ The elastic 
constants vrere talsito as preassigned valuea, i.e., the 
material to be used will be selected in advance from a 
set of available well characterized systerrjs » Closed form 
solutions for the natural frequencies, buckling loads and 
deflection are available only for simply supported lami- 
nated orthotropic plates. 

The optimisation procedure was used for the 
design of simply supported laminated plates only since 
the computer time required on IBM 70 ^+^ becomes very 



large for solving the partial differential equations for 
other boundary conditions by approximate methods. 

The following optisiisa'tion problems were solved 
in the present work ; 

Problem I s 
Given : 

1 . The available orientation angles 0^^ and the 

corresponding material properties ^Ti’ %iTi=' 

^^Ti ^i wei^t density for i = 1, 2, 

N. 

2. IT the number of available orientation angles. 

3. Lo\irer limits on the laminate fundamental frequency 

j 0 = "1 ) 2j ... k, 

4. k the number of sub-problems to be solved, 

5. Length ’a' and breadth ’b’ of the laminated plate. 

Find ; A set of t^, i = 1, 2, ..., N, such that 

N 


1 . 

Wei^t = 2 2 Pi "'"i minimum 

( 9 ^) 


i=1 


2. 

Pmn?-' Pj » j = 1, 2, k 

( 95 ) 

3 . 

ti ^ 0 for 1 = 1,2, . . . , N 

( 96 ) 


It should be recognized that equation (^) states the 
minimum weight objective, equation (95) implements the 
frequency constraint and equation (96) requires that the 
thl clones ses be non-negative. 



Problem II : 




Given : 

1 . The available orientation angles 9 ^^ and the 
corresponding material properties 

hii’ ^Lii Fi 1 = ■'> 2 

2 * Nj 5 0 = 1 , 2 , . . k, the number of available 
orientation angles . 

3 . k is the nrunber of sub -problems to be solved. 

• Lower limit on the laminate fundamental frequency, 
5. Length 'a’ and breadth ’bl of the laminated plate. 
Find: A set of t. , i = 1, 2, such that 

•J- J 

1 . VJeight =2 53 p t-. ab , — — ^ minimum (97) 

i =1 

2 - Pmn?/S . (98) 

3. t >> 0 for 1 = 1 , 2 , N (99) 

X J 

The equations (97) to (99) have the same meaning as 
equations (9^) to ( 96 ) of problem I. 

Problem III : 

Given : 

1 . The available orientation angles 9^ and the 
corresponding material properties 

G^Ti? ^LTi Pi i = 1} 2, N. 



2 . 


N the number of available orientation angles, 
3. Upper limit on the weight of the plate 

} 3 ~ "Ij 2j h 

U-. k is the number of sub-problems to be solved. 


5. Lengtii 'a' and breadth 'b' of the laminated plate, 
Find: A set of t^^, i = 1, 2, IT, such that 


1 . 

Frequency 

= maxijnmi 

(100) 

2. 

wei^it fC 


(101) 

3. 

ti ^ 0 

for i = 1 , 2, . . ., N 

(102) 


It should be recognised that equation (100) states the 
maximum frequency objective, equation (101) implements 
the v/eight constraint and equation -(102) requires that the 
thicknesses- be neii' negative . 

Problem IV : 

Given : 

1 . The available orientation angles and the 

corresponding mateilal properties ^Ti’ 

^LTi’ ^LTi <^i ^ ~ • 

2. N. the number of available orientation angles 

3 

for j = 1 , 2, . . , , k, 

3. k is the number of sub-problems to be solved. 

4-, Lower limits on the laminate fundaiiBntal 

frequency ^ • , and the critical bucldLing load c(^ . 
• 3 
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5^. Length 'a’ and. breadth ’h’ of the laminated 
plate . 

Find; A set of t . , i = 1 , 2, ...,11. 

3 


1 . 

weight = 2 

tj^ab ' . — minimum 

i=1 

(103) 

2. 

^mn h=^^'j 


( 1 o4-) 

3. 

(N ) . 

X or j 


(105) 

^h. 

t^ 0 for 

1 = 1, 2, ..., Ih 

(106) 


It should he recognized that equation (103) states the 
minimum weight objective, equation (10^) implements the 
frequency constraint and equation ( 105 ) implements the 
buckling load constraint and equation (106) requires that 
the thiclmesses be non-^negatlve . 


Problem V : 

Given 

1 . The available orientation angle 0^^ and the 
corresponding material properties 

^LTi’ ^LTi /i ^ ~ ***’ 

2. N. the number of available orientation angles 

0 * 

forj~1,2j..»jk. 

3. k is the number of sub -problems to be solved. 






4-if 

4. 

Lomov limits 

on the lam.inate fundamental 

frequency, 


, and the 

critical buckling load, o< . 


5. 

Upper limit 

on the transverse deflection 

of the 


plate, 



6. 

Length 'a' and breadth 'b' of the laminated plate. 

Find A 

set of t^, i 

= 1 , 2, . . . , 11 . , such that 

J 


1 . 

V/eight = 2 

”3 

Pi "^i minimum 

i=1 ^ 

( 107 ) 

2 . 

^mn ^ 


( 108 ) 

3 . 

(N ) > <> 

X cr"^ 


( 109 ) 

if. 

^ max 


( 110 ) 

5 . 

t^ ^ 0 for ; 

1 = 1 , 2 , . . . , H j 

(ill) 


Here' eq'uation (107) states the rainim'um -weight objective, 
equation (108) implements the frequency constraint, 
equation (IO 9 ) implements critical bucltLing load constraint, 
equation (IIO) implements, the deflection constraint, and 
equation (ill) requires that the thicknesses be non-negative. 

It can be noticed that all these five problems 
are similar to the one stated in equation (82). 





3 .2 SO LUTIOM TECT I NIQUE 

All exaiiiination of the frequency, huclfLing load 
and the deflection equations reveals that the quantities 
are expressed in terns of the material stiffnesses A. 
coupling matrix elements B. and flexural rigidities D 
etc., uiiich in turn are functions of the thickness, orien- 
tation and properties of each ply. Since B.. and D. . 

3-0 i 0 

are not linear functions of tlie ply thicknesses, these 
quantities uilO. he non-linear functions of the design 
vector i.e., ply thicknesses. The non-linear optimization 
prohlcrns are solved using the interior penalty function 
metliod in this uork. 

In tlii.s approach, augment the objective func- 
tion with a penalty terra which is small at points away 
from the constraints in the feasible region, but wliich 
’blows up* as the constraints are approached. It can be 
expressed as 

P 

0(X, r) = ]?(X)-r Z l/gj ( 2 ) (112) 

3=1 

where F is to be minimized over all X satisfying 
g- ( X ) 0, j = 1, 2, .. ., p. Note that if r is positive 

the effect is to add a positive penalty to F ( X ). This 
is because at an interior point all the terms in the sum 
are negative.’ As a boundary is approached, some g^. will 
approach zero and the penalty will "explode". It has been 



proved that the constrained minim-uun of F can he obtained 
hy finding the -unconstrained minimum of 0 for decreasing 
sequence of r *■ values, 

3.2.1 Algorithm 

The solution procedure can he described by the 
f oil owi ng steps: 

1 . Given a starting point satisfying all g^ { X ) ^ 
and. an initial value of r, minimize 0 to obtain 

2. Check for convergence of to the optimum. 

3. If the convergence criterion is not '.satisfied replace 
r by r rc -v^here c < 1 . 

4 . Take Xj^ as the new starting point and repeat 
from step 1 . 

There are a number of points to consider in apply- 
ing the method: 

a) The starting design required by step 1 is 

usually available in engineering problems, but 
sometimes finding such a point becomes difficult. 

b) A proper initial value for r m-ust be selected, 

c) The possibilities for the convergence criteria 
of step 2 are n-umerous and there are choices to 
be made. 
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d) Because of the sequential nature of the process, 
it is possible to improve the starting point for 
the third and ‘subsequent minimizations, 

3.2.2 Starting Point 

In many engineering situations, particularly in 
the structural and mechanical design areas, it is easy to 
find a point satisfying g. ( X ) ^ 0 at the expense of 

tJ 

large values of F, For example, in structural design, 
if cost or 'vrei^t of the structure is ignored, it is usually 
easy to propose many designs which fulfill the basic 
requirements of strength and rigidity for the particular 
application. In other design situations, however, the 
acceptable designs may not be at all obvious. In these 
si-tuations the initial acceptable design required by the 
interior penalty function method can be obtained as follovrs . 

Suppose an engineering assessment of the situation 
has produced tlie design X^ idiich satisfies g^ ( X^) < 0, 

2 = -1, 2, s, but which has g^ ( X^ ) ^ 0, j = s +1, 

s + 2, p. Take the k for which ( X^ ) is a 

maxim'um out of g- ( X ) for 3=s + 1,s+2, 
and temporarily define it as the objective function for 


the following problem : 



Find X such that 


gjj ( X ) -hD> miniimm 

gj ( X ) :^ 0 3 = I5 2, s 

g^- ( X ) - g^. ( X^ 0) 0=3 + 158 + 2j ... 

h*"l5h + l5 ...^n 

Whenever, during the process of solving this problem by the 

penalty function method, the value of gj^ ( X ^ below 

zero, the procedure is halted. 'The point so obtained then 

satisfies at least one more constraint than did the original 
The procedure can be repeated untill all the cons- 
traints have been satisfied and an 3?^ is obtained for 
which ( X^ ) < 0, g = 1, 2, ..., p. 

3.2.3 Initial Value for r 

The matter of selecting an initial value for the 
penalty parameter r has been discussed in the literature, 
but while some theory is available, the task is still an 
art. If r is large, the Function is easy to minimize, but 
the minimum may lie far from the desired solution to the 
original constrained minimization problem. On the other 
hand, if r is small, the function will be hard to mini- 
mize. A general rule is that if is a conservative 
design, pick so that Vg^ ^^o^ approximately 

equals F In practice, this approach usually 

yields a reasonable imtial value for r. 



3.2.4 Convergence Criteria 


4 


As the ^ - fimction is minimized for various 
decreasing values of r, the sequence of minima 

^ **“’ should converge to the solution 

of the constrained minimization problem. A means is 
needed to ascertain this convergence mthout an unnece- 
ssarily large number of minimizations. One simple crite- 
rion is to compute the relative difference 


TP' , \ ^ ''V 

Tniri^'‘^l-r ^ “^min 

__ ^ p _ ^ 

min 1 


(rj_) 


(113) 


and stop when this value drops below a certain fraction. 


where 


Another convergence test can be stated as 

1aU£ 


3 . 2.5 Unconstrained liinimization 

For the unconstrained minimization of the penalty 
function, the Variable Matric method (also known as the 
Davidon - Fletcher - Powel method) has been used. This 
method is a first order method wMch replaces the local 
hessian "by an approximate metric K^, The method of 

computing this metric completely eleminates the need 
for evaJ-Uating the second derivatives and performing 
matrix inversions, and yet the sequence of iterations 
converges quadratically to the minimum point 
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Furthermore, it can be proved that the matric which 
is improved at each iteration, converges to 

The iterative process can be stated as follov/s : 

1 • Start with an initial and an initial positive 
definite S3niiiiietric matrix (for example, the 
identit 3 ^ matrix), and set <7 , 

2. Compute X . = 3^ + .;o< s , where o(* minimizes 

4T I q q q. q 

^ check for convergence. 

3* Compute where, defining 


4. 



and 



Gq (X^^p -7F (X^) 


S S 

a 


s I 
q 


T 


q 



(H ir ) (H I ) 
q q q a 

„ T . _ 

Y H Y 


_ q q q 

Compute and repeat from 

step 2. 



51 


3,3 MJMERICAL RESULTS 


The numerical results of optimizations are 
presented in this section. The nijmerical data assumed for 
each ply in all the problems ims 

= 31.0 X 10^ psi; = 5.4 x 'l0^ psi; 

^LT ~ 0.75 X 10^ psii ^'"^LT ~ 0,28 and 
f - ,0001856 Ib/in^, 

The orientation angles assumed were 0° and 90® alternatively 
with symmetry about the middle plane. 


3 .3*1 liUminated Beams 

The common data for the two problems is 
1 = 10,0 in., X = 48,000 psi, Y :=■ 30,000 psi, S = 8,000 psi. 

Problem I : 


The number of design variables considered is R = 2 
only and the strength constraints were cubic functions of 
t^ and t2 so the problem was solved graphically. 

Objective : minimization of volume (in,^) 

Design variables : thicknesses of plies (in.) 


The values of t2 for preassigned values of t^ 
satisfying the various strength constraints are listed 
in table 3.1 »■ and the equations are shown plotted in 


figure 6. 


I.I.T. K/.NPUR 
CENTK /,L I IBRARY 

Acc. Kfo. A 
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TABLE 3 

.1 Values of t^ 

for different t^ 

b 

(inches) 

5 

Value or tp from constraint on (Inc/i 

a Normal stress 



cr SI JNormal' stress 

O' 2" Shear stress 

7 5 <^xy 

0.0 

0.098^2 

0.124-70 

0.21307 

0.005 

0.01890 

0.02992- 

0.08039 

0.010 

0.00804- 

0.014-38 

0.0+4-22 

0.015 

0.00227 

0.00719 

0.02890 

0.020 

- 0.00278 

0.00183 

0.01979 

0.025 

(-)ve values 

- 0.00320 

0.01310 

0.030 

II 

('-)ve values 

0.0074-9 

0.035 

n 

ti 

0.00235 

O.C^O 

rr 

f! 

- 0.00265 


Lines representing constant volume are also 
plotted in figure 6j and the optimum thicknesses are found 
to be 


t^ = 0.0375 in., t2 = 0.0 in and 

•3 

volume = 0.75 in. 

It can be seen that the shear stress constraint 
is active at the optimum point. This can be used as the 
optimality criterion for tliis class of beams. 
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ProlDleni II 

Though the number of design variables is I = 2 in 
this case also, it is not possible to express the deflection 
constraint in explicit form in terms of the thiclmesses 
t^ and t^, as it involved the inverse of [a] . Hence the 

problem is solved on IBM 7C^-p4- computer using the penalty 
function technique. In all, five values of r uere con- 
sidered for the minimization of - function I’Ath 10 itera- 
tions for each r (in the Davidon - Fletcher - Pouell 
method) . Three refits were permitted during each cubic 
interpolation stage. 

2 

Load acting on the beam Q = 1 Ib/in 

Upper limit on deflection, = 0,1 in. 

Objective ; Miniinization of volume (in.-^) 

Design variables : Thickness of plies (inch.) 


At starting design point ; 

X = 

0 

o 

• 

o 

• 

o 

^1 ^ 

0.07 


) = 15.3838!^ 

F ( Xo ) 

= 7.7 

V ( ) 

= 0.07^5Wf inch. 


= Wl .87559 iVin,^ 


= 123.72517 Ib/in.^ 





^xy ^ = 0-2“<-09775 X lO"® lb/in.2 

At final design point ; 

Xopt = 0.321125, O.CiK 69 V 5 -i -2 
= 0.7 X lO"^ 



= 7.361556 

^ ^ -opt > 

= 7.361383 

^ ^ ^opt > 

= 0 . 0995+5503 inch 

^ ^opt > 

= 484.32 Ib/in.^ 

XI 

0 

= 135.6 Ib/in^ 


^xy ^ ^opt ^ = 0.2^09775 X 10 "^ It/in.^ 

Throughout the optimization process, the deflec- 
tion constraint has been active. 

A comparison of the tvo beam problems indicates 
that the deflection is a critical parameter in the design 
of beams as one -tenth of the load resifLted in about nine 
times increase in volume of the beam when deflection is 
also considered as a constraint. 
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3,3.2 Laminated Plates 

The restilts for all the prohlems were obtained on 
IBM 70 ^ computer. The numerical results of optimization 
are given in this section. In all the problems, four 
values of r vrere considered for the minimisation of ja ~ 
function with 10 iterations for each' r (in the Davidon - 
Fletcher - Powell method). Three refits were permitted 
during each cubic interpolation stage. 

Problem I 

Total number of plies = 8 
Number of design variables, N = if 
Objective : minimization of weight (lbs.) 

Design variables : thicknesses of plies (inch) 

(a) iS ^ = 100.0 cycles/sec. 

At starting design point : 

= 0 . 012 , 0 . 012 , 0 . 012 , 0.012 

r^ = 1 x 10"^ 

ja ( , r^) = 0.2657^11 X 10"^ 

F ( Xq) = 0.1990018 X 10"^ 

P1I ( 


= 237.61335 cycles/sec. 
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At final design point : 

x^pt = 0,0050566, o.oo5o8c^9, 0.005054^^1, 0.00507321 

= 1 X 10“^ 

0 C = 0,Bh2(hk2 X 10“-^ 

F ( I^p^) = 0.8 !-i-01VB6 X 10“3 

P 1 I ^ V ^ = 100-3159^8 

CT 3 > f^2 = 150.0 cycles/sec. 

At starting design point ; 

X^ = 0.012, 0.012, 0.012, 0.012 

TrJ = 1 X 10“*^ 

0 ( X^, r^) = 0.2658396 x 10“^ 

F ( Xq ) = 0.1990018 X 10”^ 

p^^ ( Xq ) = 237.61335 cycles/sec. 


At final design point i 


X ^ 
opt 


= 0.00733182+, 0.00772892, 0.00765082, 0.00763795 

= 1 X 10"^ 


^ ^ %t’ = 0.1259935 X 10 

F ( X 4 . ) = 0.1258252 X 10 

^ opt 


o ( X ^ ) = 150.238555 cycles/sec. 

■2^11 ^ opt 
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(c) ~ 200.0 cycles/sec. 

At starting design point : 


= 0.012, 0.012, 0.012, 0.012 

= 1 X 10 

0 ( X^, r^) = 0.2662002 x lO^^ 

p ( ) = 0.1990018 X 10"^ 

( X^ ) = 237.61335 cycles/sec. 

At final design point 

X = o.oiocM+1, 0.010350^, O.OIOOWj 

opt . ' 

II, = 1 X 10-9 

^ ^ ^opt» V = 0*1679688 X lO"^ 

p ( X ) = 0.1678405 X 10-^ 

opt 

„ ( ) = 200.405934 cycles/s 90. 

P 1 I ^ opt^ 

(d) ^ = 250.0 cycles/sec. 


At starting design point . 

X = o.oi 5 j o.oi 5 j o.oi 5 j 0.015 

0 

= 1 X 10-^ 

1 


0 ( 1*1 ^ 
F ( ) 

P1I ^o ^ 


- 0.3026173 X 10 

= 0.2487522 X IQ-^ 

_ 297.01669 cycles/sec. 


0.0100446 
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Tile relation between minimum weight and j5 . is 
shown plot bed in figure 7(a). It can be seen that the 
minimum weight varies linearly with the lower bound value 
of the frequency. The computer time taken in each caee 
is about eight minutes . It vias observed that frequency 
constraint was active in all the problems . 

Problem II 

^ = 200.0 cycles/sec. 

i 

Objective : minimization of weight (lbs.) 

Design variables : thicknesses of plies (inch) 

(a) Total number of plies = h 

Number of design variables, = 2 

At starting design point ; 

= 0.02^+, 0.024- 

■r^ = 1 X 10’"'^ 

( X^, r^) 0.2162002 x 10“^ 

■F ( X ) = 0.1980019 s: 10”^ 

p^^ (' Xq) rr 237.61335 cycles/s eb. 

At final design ■ point : 

\pt = 0.0203055, 0.0203077 

= 1 X 10-^ 

0 ( ) = 0.1684-155 X 10"^ 

F ( Xqp^ ) = 0.1683767 x 10“^ 

( Xq ) = 201.04-6202 cycles/sec. 
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(Td) Total number of plies = 8 

tober of design variables, = k 
At starting design point : 

= 0 . 012 , 0 . 012 , 0 . 012 , 0.012 

= lx 10 “^ 

0 ( X^, r^ ) = 0.2662002 x 10“^ 

P ( X^ ) = 0.1990018 X 10"^ 

p^-, ( Xq ) = 237.61335 cycles/sec. 

At final design point i 

^opt ~ 0.0100M+1, 0. 0103504, 0.0100^7, 0.0100446 

= 1 X 10“9 

0 ( iij. ) = 0.1679688 X 10"^ 

F ( X^p^ ) = 0.1678405 X 10"^ 

^11 ^ ^opt ^ ' = '200.40593^ cycles/sec. 

(c) Total number of plies =12 

Number of design variables, = 6 

At starting design point ; 

Xq = 0.008, 0.008, 0.008, 0 . 008 , 0.008, 0.008 
- lx 10“^ 

0 ( X^, r., ) = 0.3495335 X 10"2 

F ( X^ ) = 0.1990018 X 10“^ 

p^^ ( ) = 237.61335 cycles/sec. 
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At final design point : 

\pt = 0 . 00697087 , 0 . 00658763 , 0 . 00676488 , 0 . 00649378 , 
0 . 00694369 , 0 . OO667 180 ' 

\ = 1 X lO”^ ■ 

^ ^opt’ ^ = 0.1679376 X 10“2 

F ( ) = 0.1676285 X 10~2 

^11 ^ ^opt ^ ~ 200.1527876 cycles/sec. 

(d) Total number of plies = 16 

Number of design variables,!^ = 8 

At starting design point' : 

Xq = 0 . 006 , 0 . 006 , 0 . 006 , O.OO6, 0 . 006 , 0.006, 

0.006, 0.006 

r^ = 1 X 10"^ 

0 ( X^, r^ ) = 0.4662002 x 10 "^ 

F ( X^ ) = O.199OOI8 X lO”"^ 

p^^ ( X^ ) = 237.61335 cyclers ec 

At final 8osign point : 

x^p^ = 0.00507456, 0.00507839, 0.00506896, 0.00506894, 
0.00506783, 0.00507637, 0.00506783, 0.00506783 

= 1 X 10"^ 

0 ( X^p^, iij. ) = 0.1685403 X lO'’^ 

F ( X^p^ ) = 0.1682009 X 10“2 

P1I ( X^p|. ) = 200,83628 cycles/sec . 
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(e) Total n-umlDcr of plies = 24 

Nurater of closing variables^N^ =12 
At starting design point ; 

= 5”. 05 x 10 5.05 X 10“^, 5'. 05" X 10“2, 5,0^ x; 10 “*^^ 

5.0? X 10 5.05 X 10“^, 5.05 X 10"^, ^.05 x.10“^y 

5.05 X io*-3, 5.05 X 10-3 , ^ ^q- 3. 

r^ = 1x10“^ 

0 ( Xq, r^) = 0.7266872 X 10-2 
( “^0 ) = 0.2512397 X 10-2 

^^11 ^ \ ^ “ 299.9868 cycles/sec. 

At final design point i 

^opt "" 3.37050 X 10-3, 3.37050 X 10-3, 3.37060 X 10-3, 

3.37010 X 10-3, 3,37010 X 10-3, 3.37010 X 10t3, 
3.37050 X 10 - 3 , 3.37000 X 10 - 3 , 3.370^0 X 10 - 3 , 

3.37010 X 10-3, 3.37010 X 10-3, 3.37000 X 10“3. 

-'lx 10 -^ 

0 ( i|j_) = 0.1684815 X 10-2 

F ( ) = 0.1676713 X 10-2 

p^^ ( X^p^ ) = 200.203916 cycles/sec. 

The relation between minimiim weight and. the number 
of plies is shown plotted in figure 7Cb). It can be observed 
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that the change in number of plies has a negligible effect 
on the rainimura weight for a preassigned lower bound on the 
natural frequency. The computer time used for the above 
problems varied from approximately seven minutes to thirty 
five minutes depending upon the number, of design variables , 
The frequency constraint was active in all the problems. 


Problem III 

Total number of plies £= 8 
Number of design variables , N = 4 

Ob,icctive •* maximization of frequency’ {eycles/sec. ) 

Design variables : thicknesses, of; plied-' Gin. ) 


(a) W..J - 0.8 X lO"^ lb. 

At starting design point : 

X = ^.0 X 10“3, 4.0 X 4.0 X 10"^, 4.0 x 10~^ 

0 

. r.^ = 0.3 X 10'^ 

0 ( X^, r^ ) = -29.66984 

F ( X^ ) = -59.84564 

W ( Xq ) = 0.6633392 X 10"^ ib. 


.At final desi^^n point ; 



4 .79!" 60 X lO'^, 4.78922 X 10'“-'^ j 4.79181 X lO*"-^, 
4.80043 X 10”^ 




Df 


0.3 X 10 



-71 .70992 
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^ ' 4pt- ) = 
w ^ 4pt ) 


-71 .73993 

0.79518006 X 10 “^ Id, 


(6) = 0.125 X 10“^ Lb, 

At starting design point ; 


Xq = 5.0 X 10“3, ^.0 X 10"^, 5.0 X 10"^5 5.0 X 10' 


ro 


= 0.5 X 10 


r1 


0 ( Xq, r., ) -34.65830 

F ( Xq ) = -74.80682 

W ( Xq ) = 0.08291725 X 10"^ lb. 

It final design point : 


Xqp^ = 7.47819 X 10“^, 7.49888 X lO-^, 7.48639 x 10”-', 

7.47108 X 10"^ 

% = 0,5 X 10“^^ 

( X^p'e 5 1-4 ^ =-11 1-9317 

F ( X^p^ ) = -111.9654 

W ( X^p^ ) = 0.12410482 X 10“^ lb. 

(c) = 0.17x 10"^ in. 

At starting design point ; 

X^ = 7.0 X 10”^, 7.0 X 10“^, 7.0 X 10"3, 7.0 x-10”^ 

r.^ = 0.2 X 10"^ 


0 ( X^, r^) 


-93.23791 
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(d) 


) = “10^.7295 

w C Xq ) = 0.ll608if33 X 10“^ lb. 

At final design point t 

^opt “ 1*02802 X 10 1.01796 X 10“^y 9.951O3 X lO'^ 

1.03381 X 10‘2 



= 0.2 X 

10“^ 

0 ( 

^opt» ' 

= -152.4068 

d' ( 

X , ) 
opt 

= -152.4179 

W ( 

X , ) 
opt ^ 

r= 0.1 6994308 X 10 

\ 

= 0.2 X 10" 

lb. 


At starting design point : 

Xo = 9.0 X 10 9.0 X 10“^, 9.0 X 10“^5 9.0 X 10“^ 

= 0.15 

0 i ) = -67*39446 

S' ( Xq ) = -134.6523 

^ \ ^ = 0.1493514 X 10"^ lb. 

At final design point t 

^opt ~ 1.186+3 X 10 1.20599 10“'^, 1.20009 X lO"^, 

1.18610 X lO-^ 

^ = 0.15 X lO”-^ 
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^ ^opt’ = “178.67Cif 

^ ^opt ^ = -178 .7365 

^ \pt ^ = 0.198115224- X 10-2 

(e) = 0.25 X lO’"^ jLls, 

At starting design point : 

Xq = 1.0 X 10“2^ 1.0 X 10"2^ 1.0 X 10"2, i .0 x lO"^ 

r^ = 0,2 

75 ( Xq, r^) = -62.01957 

F ( X^ ) = -1^9.6136 

W ( X^ ) = 0.16583^(-75 X 10”2 

At final design point : 

Xq = 1.4-92Mf X 10“2^ 1.4-953A x 10~2^ 1.4854-8 x 10"2, 

1.4-8712 X 10“2 

= 0.2 X lO-^ 

^ ^ ^opt’ ^ ^ ^ -222.8677 

F ( X^p^ ) = -222.9387 

^ ^opt ' = 0.24-71096 X 10“2 113. 

The relation between maximron frequency and the 
design weight is shown plotted in figure 7(c). It can be 
obserA/Bd that the maximum frequency varies linearly Ad.th 
the preassigned design weight as seen in problem I. The 
computer time used for each problem is approximately 
eight minutes , 
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^ = 200.0 cycles/sec. 

0( = 1700.00 Ib/in. 

Objective ; minimization of wei^t (lbs.) 

Design variables ; thicknesses' of plies (inch) 

(a) Total number of plies = V 

Number of design variables, = 2 

At starting design point : 

= 5.07 x 10~2^ 5.07 X 10"^ 
r^ = 0,7 X 10”^^ 

0 ( Xq, ) = 0.6308014 X 10 "'^ 

P ( ) = 0.4187329 X 10 “^ 

( X^ ) = 337.77443 cycles/sed. 

(Ny_)^^|_ = 2390.3541 Ib/in. 

i ^0 

At final design point : 

Xopt = 4.34760 X 10 "^, 4.3476O X 10 ”^ 

= 0.5 X lO"*^ 

0 ( X^p^, i^) = 0.36102722X 10“^ 

F ( ) = 0.36(4920 X 10"^ 

P1I ( ) = 327.2304 cycles/sec. 

(N ) I _ 

x^crs V 

! ""out 


= 1727.2390 Ib/in, 



(Id) Total number of plies = 8 

Number of design variables, N^, = 4- 

At starting' design point ; 

= 2.525 X 10 *"^, 2.525 X lO"^, 2.525 X lO"^, 

2.525 X 10”^ 

r^ = 0.5 X 10 "^ 

( Xq, r^ ) = 0.122i+86l-x- IQ”'’ 

f = 0 . 4-187329 X 10“^ 

P^-j ( Xo ) = 377.774-4-3 cycles/sec. 

(N I _ = 2390 . 354-1 Ib/in. ■ 

A UX £ -y 

i ^0 

At final design point : 


Xopt “ 2.1728 X 10“^, 2.17378 X 10 “^, 2.174-71 X 10 
2.17371 X 10"^ 

11 ^ = 0.5 X 10“^ 

0 ( X^pt’ ^ 0.3617100 X 10“2 

1 ( X^p-t ) = 0.3604-835 



= 1525 . 1304 - Ib/in. 
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(c) 


Total Number of plies =12 
Number of design variables, N^ = 6 


At starting design point ; 

= 1.7 z 10~^, 1.7 X 10"^, 1.7 X 10“^, 1,7 X 10"^, 

1.7 X 10“^, 1.7 X 10~^ 


7 ^ = 0.^ X 10“^ 


0 ( X^, ) 


0,6006799 X 10' 


-2 


F ( ) 


= 0 . if 228788 X 10 


-2 


( Xq ) = 381.51^78 


^Vcr 


25+62,0600 


X. 


At final design point ; 

\pt = 1*^355 X 10“^, i.4Wf8 X 10"^, 1.5+5+366 x 10~^ 
1.5+5+395 X lO"^, 1 . 1 + 545+17 X 10**^, 1.5+1+396 X 10“^ 


= 0.5 X 10“® 


0 ( X 


opt* 




( ^opt > 
^11 ^ ^opt ^ 


^^x^cr 


!^opt 


0.3995+827 X 10 
0.3591885 X 10 “^ 

325+. 05545 + 


1508.75259 
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(d) Total Nimiber of plies = 16 

N-uinber of design variables j = 8 

At starting design point : 

= 1.3 X 10~2, 1.3 X 10"^, 1.3 X 10"^, 1.3 X 10"2, 

1.3 X 10"^, 1.3 X 10"^, 1.3 X 10“^-, 1,3 X 10“^ 

= 0.5 X 10^ 

( X^,’ r^ j = 0.3520llf3 X 10"”^ 

F ( X ) = 0.^1:311705 X 10“^ 

O 

) = 388.995^5 cycles/sec. 

- = 2609.7^53 iVin. 

IX CX -y 
■^O 

At final design point : 

X^p^ = 1.0827 X 10”^, 1.08271 X 10**^, 1.08309 X 10“^, 

1.08315 X 10“^, 1.08398 X 10“^, 1.08397 X 10“^, 

1.08358 X lo"^, 1.08358 X 10“^ 

= 0.5 X 10“'^ 

0 ( X^p^, 11 ^) = 0.3637399 X 10 “^ 

^ ^ %t ^ " 0.3593131 X 10"2 

■ Pl1 ( \pt ^ ■ "" 324.167V cycles/sec. 

= 1510.3240 iVin. 

opt 
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(e) Total numter of plies = 20 

‘ Humber of design variables^H^ =10 
At starting design point : 

= 1.01 z 10"^, 1.01 X lO"^, 1.01 X 10"^, 1.01 X 10“^, 

1.01 X 10"^, 1.01 X 10'^, 1.01 X 10"^, 1.01 X 10“^, 
1.01 X 10“^, 1.01 X 10“^ 
r^ = 0,5 X 10”^ 

0 ( Xq, r^ ) = 0.5383276 X 10"'’ 

F ( X^ ) = 0.4187329 X 10"^ 

p^^ ( X^ ) = 377.77443 cycles/sec. 

- = 2390.354 Ib/in 

Ji, X I -y 

'^o 

At final design point : 

^opt 8.70907 X 10"8, 8.70907 X 10"^, 8.70854 X 10"^,. 

8.71445 X 10"8^ 8.70679 x 10"8, 8.72347 x 10"^, 
8.71662 X 10"8, 8.71035 X 10"8, 8.72122 X 10"8, 

8.72304 X 10"8 


= 0,5 X 

10 "^ 

^ ^ 4 pt’ V 

= 0.367 2401 ::-x" 10 "’^ 

^ ‘ 4 pt > 

= 0.3612819 X 10"^ 

^11 ^ ^opt^ 

= 325.9430 cycles/sec. 

fJ'xbr I, 

= 1535.2872 Ib/in. 

■^opt 
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The relationship between the minimum weight and 
the number of plies is shown plotted in figure 7(<i)» It 
can be observed that the minimnm weight varies with the 
change in number of plies for preassigned values of lower 
bounds on natural frequency, and critical buchling load. 
The optimum number of plies seems to be 1^ (tvrilce N). The 
coi^puter time used for each problem varied from about 
seven minutes to thirtj'-five minutes depending upon the 
number of design variables. The critical buckling load 
constraint was active at optimum point in all the ceases. 


Problem V 


P - 200.0 cps. 

= 1,500 Ib/in. 

"max = ^ 

Objective » minimization of weight (lb.) 

Design variables : thickness of plies (in.) 


(a) Total number of plies = 8 

Number of design variables, N^ = ^ 

At starting design point ; 

Xq = 2.5 X 10“^, 2.5 X 10“^, 2.5 X 10"^, 2.5 X 10 

r^ = 0.3 s lO”-^ 

0 ( X^, r^) = 0.473791 2; 

F ( X^ ) = 0. 4145870 X 10"^ 
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Cb) 


P11 ^ ^0 ^ 


(V=r 


X 


¥ ( X ) 

max o 


37^.03409 cycles/sec. 
2320.05^1 lb ./in. 

0.38836152 X lO*’^ in. 


At final design point j 


X 4, 
opt 


2.9O605 X 10“^, 1.99260 X 10 "^, 2.43860 X 10“^, 

2.32236 X 10 "^ 



= 0.3 X 

00 

1 

0 

^ ^ \pt> 

= 0.4005675 X 10"2 

^ <■ %t ^ 

= 0 . 40 C 4750 X 10 “^ 

P1I ( ^0 ^ 

= 361 .3024 cycles/sec 

«2or 

X 

= 1941.111 lb ./in. 


0 


w ( 

max 

^^0 ^ 

= 0,39713621 X 10"^ in 


Total Number of plies =12 
Number of design variables, N^ = 6 
At starting design point ; 


Xq = 1.7 X 10 “^, 1.7 X 10 ~^, 1.7 X 10 ”^, 1.7 X 10 “^, 
1 .7 X 10 “^, 1 .7 X 10 "^ 

r^ = 0.25 X 10“^ 

0 ( X^, r^ = 0.5267113 X 10“^ 


F ( X^) 


'■2 


0.4228788 X 10 



7 ^- 


(c) 


w ( X ) 

max 0 

= 0.3933030 X 10”^ in. 

P11 ( 

0 ) 

= 381.51478 cycles/sec. 

(Vor 

X 

0 

= 246 2.060 lb ./in. 


At final design point ; 


Y 


opt 


2.00051 X 1.3if472 X 10”^, 1.76682 X 10“2, 
I.52I63 X 10“^, 1.62367 X 10"^, I.598OO X lO"^- 




0.25 X 10' 


-8 


^ ^ %t’ = o.i+089727 X 10 


-2 


F ( X , ) 
opt 

^max ^ ^opt ^ 


= 0.4085897 X 10"^ 

= 0.3999654 X 10~^ in. 


P11 ( ^opf ^ - 368.6234 cycles/sec. 


L = 2220.819 It ./in. 

cx I ^ 

■ opt 

Total i-limber of plies = 16 


Number 


of design variables , = 8 


At starting design point ; 

= 1.3 X 10"2, 1.3 X 10“^, 1.3 X 10“^, 1.3 X 10~^, 
1.3 X 10"^, 1.3 X 10“2, 1.3 X 10"^, 1,3 X 10*“^ 

r^ = 0.15 X 10 ~^ 

/ ( X^, r^) = 0.4818080 x 10"^ 
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P (X^) = 0.4311705 X 10-2 

\ax ^ ^ “ 0.38544437 z 10”^ in. 

^1 1 ^ ^ ~ 388.99545 cycles/sec. 

y 

i 

= 2609.7453 lb ./in. 


At final design point i 

^opt = 1.69641 X 10“2^ 1.02909 X 10“2^ 1.50000 X lO^^^ 
1.17626 X 10 1.18935 X 10“^, 1-. 07873 X lO^^j 

1.12785 X 10 ”^, 1.11559 X 10 “^. 

\ = 0.15 X 10 “^ 


(d) 


^ ^ ^opt’ 

^ ^ %t ^ 

\ax ^ ^ 

P 1 I ^ %t ) 
^VcrL 

Total IT-umber of 


= 0,4110638 X 10 “^ 

= 0.4109910 X 10"2 

= 0.3977044 X 10"2 in. 

= 370.7898 cycles/sec 

= 2260.2045 Id ./in. 

plies =20 ^ 


ITmber of design variables, = 10 
At starting design point s 


Xq - 1.1 X 10"2, 1.1 X lO’*^, 1.1 X 10“^, 1.1 X 10”^, 

1.1 X 10'^, 1.1 X 10“^, 1.1 X 10 “^, 1,1 X 10“^, 
1 .1 X lO"^, 1 .1 X 10 “^ 

= 0.25 X 10 "^ 
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0 ( Xo> 1-1 ) 

= 0.6856721 X 10“^ 

IM Xo ) 

= 0,4-56874-9 X 10**^ 

iiax ^ ^ 

= 0.331724-69 X 10"2 in 

P11 ^0 ^ 

= 4-12.18558 cycles/sec. 

1 

( ) 1 
x^cr ^ 

i 

= 3 104- ,8663 lb, /in. 


At final design point ; 

^opt “ 1*30500 X 10 7.56153 X 10'’-^, 1.16610 X 10*“^, 

8 . 4-^384 X 10“^, 1.074-92 X 10“^, 9.11366 X lO*"^, 

1.01533 X 10“2^ 9.56735 X lO''^, 9.80984- X 10"3j 

9 . 74-539 X 10“^ 

11^ = 0,25 X 10 “^ • 

^ ^ ^opt’ 0.4-1if294-9 X 10“^ 

^ ( ^opt ^ ~ 0.4-139862 X 10 ^ 

\ax ^ ^opt ^ = 0.3980690 X 10"2 in 

P-11 ^ ^^opt ^ “ 373.4-92 cycles/sec. 

^ T\T "S 

>- cr y- = 2309.982 lb, /in. 


(N ) 

X cr 


The relationship bet^jeen the miniraiUK veight and 
the number of plies is slioini plotted in figure 7(e). It 
can be observed that the minimum weight varies with the 
change in number of plies for preassigne^i vaAues of lower 
bounds on natural frenuenc;^ and critical buchling load and 
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fclie upoor "bound on transvors© deflsction. The optlinmii nuiiiber 
of plies seems to he 12 . The computer time used for each 
prold.ori varied from hine'. -.minutes to thirtjrfive minutes 
derH).ocLlnc upon the number of design variables. The trans- 
verse deflection constraint ms active at optimum point 
in a3.1 the cases , 





CHAPTER 'IV 


EXPERIIviENTAL STUDY 

The aim of this work is to find the natural 
frequencies of vihration of fibre reinforced orthotropic 
laminates at different fibre orientation angles because 
closed forra solutions are not aTaila.ble for them. 


h . 1 EXPERIMENTAL CTiARACTERIZAT IOH 
OF A LAMEN A 

In order to find the natural frequencies of com- 
posite plates theoretically, five intrinsic macroscopic 
properties E^, E^, "^TL ^LT ^ basic lamina 

or ply are required. 

In theory, discrepancies e^cLst between the various 
micro -mechanical expressions for predicting the fibre 
elastic constants. Therefore, it is necessary, certainly 
for the more important structures, that these quantities 
be based on a phenomenological approach and determined 
experimentally. A simple teclinique for accomplishing this 
is discussed in tliis section. 

At lee.st three test are required to deteririne 
the desired elastic properties in oxrj direction Q from 
the following theoretical information ^ 
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E. 


4 - 

cos 0 


+ ( 




“ Q 9 + ^ — ~ (115) 


E, 


E, 


. ‘If 

sin 0 


E, 


T 


+ ( 


sin^ @„gQg2 Q ^ C ^ s , 9 d-i^) 

^LT ®T 


1 


12 


cos^ 20 


G. 


+ ( 


LT 


1 -i^T m 1 + Vrp^ ^ 

—M + — sin^ 2 9 

El E^ 


( 117 ) 




12 


E 


V, 


^ ^ 1 ( ^LT ^ '^^LT 1 ^ - 2 Q 

E-l If'-E “ n J sin 2 9 


21 


E„ 


G 


LT 


( 118 ) 

where the suhscripts L and T denote the longitudinal and 
transverse directions respectively. 

These tests can be simple tensile testsas shown in 
figure .1C . 

The test of figure 10(a) can be used to determine 
the elastic constants E^ and Ey measuring the 

load P, the cross-sectional area A, and the strains 
and £ 


E 


^1 


1 “ - A P 


and 


■1 


^ 1 


12 


■1 


( 119 ) 


can be established, since 0 = 0*^ for this case, equations 
(1l5) to (116) give, 


1 

E. 


1 


1, 


or 


E. 


E. 


A e 


( 120 ) 
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and 




^ ) 
or 


LT 




12 




( 121 ) 


Since is a tensile strain (positive) and 6^ a 

compressive strain (negative), is positive in sign 

in equation (121). 

In a siniilar manner, the test of figure 10(1) 
where 9 =90°, yields, 


E 


T 




and 





( 122 ) 

(123) 


A check on the accuracy of the measurements can he made at 

this point, since the equality, Ej^ = E^ ^LT’ 

hold. 


The final elastic constant 0 ,^ can he determined 


hy the test of figure 10(c). 
that 0 = ^5° and 

TD 


LT 

If the test is made such 


E 


1 


A e 


1 


is computed, the value of can he determined from 

equation (115) as ; 


1 


E. 


i.e. 




If 


2 


E, 


+ 


LT 


E. 


^LT 


=L 


2 V 


G 


LT 


E. 


E, 


Ir 


+ 


LT 


E. 



01 


if .2 EXPERII^lELr rA L SET-UP FOR ^EA SlTBTNrT NATUEAL PREODlTiClES 
OF A LAl-ilMT ED PLAIE ^ " 

The natural frequencies of vibration for simply 

supported square plates were determined experimentally. 

To simulate supported boundaiy conditions, initially balls 

of 1/4" diameter -vjere used. These balls were placed at a 
1 

distance of 1^ along the edge on both sides of the plate. 
The lov/er set of balls was placed in holes in an angle 
structure fixed to the ground and the upper set of balls was 
held in the holes of m.s. flats with the help of grease. 

To varify the boundary conditions simulated, the natural 
frequencies of an aluminium plate were found on the test 
set~up. A tightening pressure of about 1/2 lb. -in., f/iiii 
the help of a torque wrench, was applied for the aluminium 
plate and the results obtained were found not to be matcliing 
with the theoretically calculated values , Also the repro~ 
ducibility of the results could not be achieved. 

Hence another attempt was made by replacing the 
balls by a 0.2 in. diameter bright bar, as sho\m in 
figure 11 . The tightening pressure applied was just 
sufficient to make the rod just touch the test plate. Tliis 
was verified by making sure that a razor blade does not 
slip in between the rod and the plate on either side. 

With this arrangement, an excellent agreement 'between the 
experimental and the theoretical values, vri.th a very good 
reproducibility of the results for aluminium plate, was 
obtained. 
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The test plates were excited indirectly mth 
acconstic excitation from a loud speaker. A sinusoidal 
input from a RC oscillator through a power araplifier was 
fed to the loud-speaker which imparted the excitation 
energy to the plate throng sinusoidal sound-x/aves . An 
electronic counter \/as used to note the frequency of the 
input signal. The block-diagram representation of the 
set-up is given in figure 10. ■ . 

A small piece of magnetic material (razor blade) 
was fixed to the composite plate with the help of cellulose 
tape and a contactless pick-up was used to find the fre- 
quency of vibration of the plate xdiich was found to be 
matching \d.th the input signal at resonance. 

The mode shapes were detected visually by spray- 
ing very fine sand on the plate surface . The final 
adjustment of the input frequency was made by displaying 
the signal of the contactless pick-up on an oscilloscope 
to ensure that the signal has a maximum amplitude. 


If. 3 TEST SPECimiT 

The vibration tests (for natural frequencies) 
were conducted on uni -and bi-directionally reinforced 
composite plates. In bidirectional plates, five layers 
of the E-glass fabric and a polyester laminate cured with 
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MEK peroxide and cobalt octate were used but the orienta- 
tion of all of them was same with respect to the plate 
axis. The orientations tested \rere 0°, 30 ° ancL 

since the angles between ^5° and 90 ° v/ill be just the 
repetition of these only. The plates were terminated to 
a size of 9»2'' x 9 •2" and the thickness was approximately 
0.034-". 

Three unidire ctionally reinforced plates irere 
fabricated from E-glass roving and Cl 230 epoxy resin 
(cured with lO^o HY 951 hardener). These plates were 
almost identical in thickness and in volume concentration 
of reinforcement. They were cut to the same size and 
glued together using the same resin as the bonding agent; 
the 90 ° plate was kept in the middle layer. The cured 
plate was trimmed to 9 *2" x 9*2" size and tlie average 
thickness of the plate was 0.1257''- The plate was so 
thick that only the fundamental frequency could be excited. 

k .4- EXPERIIIENTAL BESULTS 

In this section the results obtained experi- 
mentally are listed. Since the fibre reinforced compo- 
sites are not readily available commercially, a limited 
number of tests were performed. 



4.4.1 Elastic Constants 

Hie tests as described iii article (4.1) were 
carried out on an Instron Testing tlacliine. The strains 
were measured with the help of strain gauges mounted on 
the test specimen. The load vs strain is shoim plotted 
in figures 8 and 8(a) for uni - and bi-directionally 
reinforced specimens respectively and the elastic cons- 
tants found are given below. 

Eor unidirectional specimen : 

= 2.137 X 10^ psi; E^ = O.3OI x 10^ psi; 

= 0.257 X 10^ psi; = 0.28 and 

f = 0,00439 lb /in? 

For bidirectional specimen : 

Ej^ = 15.27 X 10*^ psi; E^ = 14.5 x 10^ psi; 

= 0.189 X 10^ ; = 0.1475 and 

P = 0.002083 Ib/in.^ 

4.4.2 Natural Frequencies 

The natural frequencies of vibration of bi- 
directionally reinforced composite plates at different 
orientation- angles are given in Tables (4.1) to (4.4), 
The size of all the platos lathin the simply supported 
edges was 9.015" x 9,017" and the overall plate dimen- 
sions were approximately 9«2" x 9 •2", 



For 0° plate : Thickness of the plate = O.O371 inches 

TAI.'LE 4.1 Natural Frequencies (cycles/sec.) 



4 729.0 


The mode shapes of vibration are shown in Figure 13. 

For 15 "° plate : Thickness of the plate = 0,0343 inch 

TABLE 4,2 Natural Frequencies (cycles/sec.) 




T 


5 

n!^<^ 

t 


0 

2 



i) 


0 

0 

56.4 


246.7 

635.2 

2 

261 .5 


457.^ 

- 

4 

670.3 


- 

- 


The mode shapes of vibration are shown in Figure 14 . 
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For 30 ° plate : Thickness of the plate = 0.0341 inch 
TABLE 4.3 'Natural frequencies (cycles/sec .) 


m — ^5 


T 


1 

n 0 

^ \ 0 

0 


2 

5 4 


D 


3 

0 

59.5 


235.0 

593.8 

2 

250.9 


463 .4 

- 

4 

615.3 


- 



The mode shapes of vibration are shown in Figure 15» 
For 45 ° plate : Thickness of the plate = 0.0342 inch. 

TABLE 4.4 Natural frequencies (cycles/sec.) 


nP< 3 

r 

0 0 

0 

1 

1 

5 2 

3 

0 

62.8 

- 

249.2 

1 

- 

224.0 

- 

2 

258.0 

- 

280.0 


The mode '.shapes of vibration are shown in Figure 16. 

The unidirectionally reinforced laminated plate 
was very thick and the excitation energy required for 
getting higher modes of vibration was much more than the 
energy of the available source of excitation. 
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The fundamental frequency of the plate was found 
to he 157*5 cycles/sec. 

Theoretical fundamental frequency of the plate 
for varying ply thickness ijith tMckness of all three 
plies same, is shoim plotted in Figure 17. 

h . 5 DlSgjSSlOU OF E.TPERIi.IEiTTAL HESTJL TS 
¥,5.1 Nodal Patterns of Orthotropic Plates 

The nodal patterns of the laminated plates -which 
behave as orthotropic plates may be divided into three 
classes; -one consisting of parallel nodal lines, one 
involving non-parallel nodal lines and one having compoun- 
ded normal modes. 

For the simply supported boundar3^ conditions, 
the mode shapes are comprised only of the rectangular 
antinoda .1 regions with nodal lines running parallel to 
the plate edges . The parallel nodal lines attribute to 
the existence of (exact) mode shapes for which a separa- 
tion into a sine wave in tlie x - direction and a irnve- 
forra in the y - direction indeijendent of x is possible. 

In the steady -state vibration considered, 
certain theoretical modes were not detected. These modes 
have an odd n-umber of nodal lines in the x - and y - 
directions. This phenomenon occurs when two opposite 
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edges have similar edge conditions and the plate cannot 
vibrato in the modes having odd number of nodal lines 
parallel to these edges under the imposed displacement 
excitation. This also can be explained physically by the 
method of excitation, and the symmeti^^ of a plate geometry, 
elastic properties and boundaiy conditions about the line 
bisecting the distance between these two edges. The 
equivalent uniform external force induced by the applied 
imposed displacement is also symmetric. Since the mode 
shapes depend exclusively on these factors , the plate can 
vibrate only in the symmetric mode shapes \diich have only 
even number of nodal lines , 

^♦7.2 Nodal Patterns of Anisotropic Plates 

No attempt was made to predict the nodal patterns 
of anisotropic plates. The effect of anisotropy 

causes an inclination of straight, otherwise parallel, 
lines or a skemng and bending of patterns. This was 
predominately clear for ^5° orientation. 



CHAPTER V 


CONCLUSIOHS AI'ID RECOLHEHDATIONS 
5.1 COHGLUSIOHS 

1 , Tlie present study shows the feasibility of 
obtaining theoretical optimLmi parameters of 
laminated composite beams and plates subjected 
to multiple behaviour constraints. 

2. The curves plotted in figure 7 for the five 
theoretical optimization problems show that the 
frequency of vibration of laminated orthotropic 
plates varies linearly with a change in the 
v/eight of the plate i,e. the thickness of the 
plate, and is fairly insensitive to the change 
in number of plies. Hovrever, the critical 
buciaing load and the transverse deflection 
depend upon the number of plies . The minimum 
weight of the plate for a constant bucl-d-ing 
load decrease first and then increases \d.th an 
increase in the number of plies, A similar 
trend is observed in the case of bound on the 
transverse deflection. 

3. The experimental study indicates that the 
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square ortliotropic plate Increases as the fibre 
orientation angle varies from 0° to ^5°. ilo 
definite pattern has been observed for the 
higher natural frequencies, 

5.2 RECOIlflElDATIOIIS ^OR FtJTUBE WORIC 

1 . The optiiirLzation of laminated composite plates 
can be attempted by talcing the fiTn?©- -orientation 
angles of the lamina as design variables . This 
requires the development of analysis procedures 
for predicting the natural frequencies and 
critical buclcling loads of such plates . 

2. The behaviour of the critical buckling load and 
the transverse deflection can be predicted 
theoretically for constant number of plies and 
varying values of the weight of the plate, 

3. The optimization of composite shell structures 

can be done in a similar manner. ’ 

It is recommended to employ another means of 
excitation so that all modes and various . 
cornpoimded nodal patterns consisting of same 
constituent modes can be detected. This can 
be achieved through the use of point exciters, 
which are capable of exciting the plate at any 
location. 
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5’. In order to verify the reduced effective 

flexural stiffnesses, a further experimental 
study should he made hy testing a series of 
unbalanced laminates and by raeasuring the 
plate surface strains or stresses. 
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Rg. I. Ply Orientation Relative to Structure ,Axes 



Fig. 2 Laminated Plate Cross-Section 
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Fig 5 Laminated Element Subjected To Moments And Forces 
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Fig. 8 Load- Strain CurvesTbr The 'pnidirectional 







Fig. 10 Block Diagram Of The Experimental Set tip 
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Fig. II Cross-Section Of Test Fixture 
( A Schematic Representotion) 
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m/n = 4/0 



Fig. 13 Mode Shapes For 8 = 0** Plate 
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Fig. 16 Mode Shapes For 9 = 45^ Plate 









Fig. 17 Theoretical Variation Of Frequency With 
Thickness, 
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Fig. 18 Photograoh Of The FynArlmon+rti 


